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Ferromagnetism in tlie t-t' Hubbard model is investigated on a square lattice. Correlation effects 
in the form of self-energy and vertex corrections are systematically incorporated within a spin- 
rotationally-symmetric scheme which explicitly preserves the Goldstone mode and is therefore in 
accord with the Mermin- Wagner theorem. Interplay of band dispersion and correlation effects 
on ferromagnetic-state stability are highlighted with respect to both long- and short-wavelength 
fluctuations, which are shown to have substantially different behaviour. Our approach provides a 
novel understanding of the enhancement of ferromagnetism near van Hove filling for t' ~ 0.5 in 
terms of strongly suppressed saddle-point contribution to the destabilizing exchange part of spin 
stiffness. Finite-temperature electron spin dynamics is investigated directly in terms of spectral- 
weight transfer across the Fermi energy due to electron-magnon coupling. Relevant in the context of 
recent magnetization measurements on ultrathin films, the role of strong thermal spin fluctuations in 
low dimensions is highlighted, in the anisotropy-stabilized ordered state, by determining the thermal 
decay of magnetization and Tc within a renormalized spin-fluctuation theory. 



I. INTRODUCTION 

Recent magnetization measurements in ultrathin Fe 
and FeCo filmsii^ have highlighted the role of dimen- 
sionality and thickness-dependent magnetic anisotropy in 
controlling thermal excitation of spin waves and thereby 
determining the temperature fall off of magnetization and 
hence the transition temperature Tc- The temperature 
dependence of spontaneous magnetization was found to 
be well described by Bloch's law, and with decreasing 
film thickness the spin-wave parameter was found to in- 
crease significantly compared to the bulk. This stabiliza- 
tion of ferromagnetic (FM) order against thermal fluc- 
tuations by magnetic anisotropy in ultrathin transition- 
metal films, observed even for a single monolayer of Fe 
on different substrates (Au,W))^i^ is relevant not only 
from the practical viewpoint of suitability of ferromag- 
netic materials for magnetic data storage applications^'' 
but also in the context of the Mermin- Wagner theorem 
which implies vanishing for two-dimensional isotropic 
magnets due to divergent contribution of long- wavelength 
spin fluctuations and finite but low Tc for small magnetic 
anisotropy. 

With regard to the role of anisotropy on the thermo- 
dynamic properties of ultrathin films, most of the earlier 
theoretical investigations have been carried out either 
within the Heisenberg modelr'^^^'^^ where the band- 
ferromagnetic nature of transition metals is ignored, or 
within the Hubbard model^'^ however, without taking 
into account the collective excitations which are of crucial 
importance in low-dimensional systems. 

Theoretical investigations of collective excitations 
in the context of transition-metal ultrathin films 
have been carried out extensively in the FM ground 
state of the Hubbard model including realistic band 
structureji^'i^ii^'i^ However, transverse spin fluctuations 
were studied in the random phase approximation (RPA), 
which neglects correlation effects and overestimates the 



stability of the ferromagnetic state. 

Due to the presence of various competing ground 
states in the t-t' Hubbard model, the issue of ferromag- 
netic ground state on a square lattice is itself a prob- 
lem of considerable recent interest. The special favor 
of the square lattice for the antiferromagnetic ground 
state at van Hove filling due to the nested Fermi surface 
is suppressed in the presence of next-nearest-neighbor 
(NNN) hopping (t'), which paves the way for various 
other instabilities including the ferromagnetic ground 
state in a narrow density range around the van Hove 
filling if t' is sufficiently large, as reported recently us- 
ing a variety of approaches. These include the Hartree- 
Fock (HF) approximation,^^ quantum Monte Carlo 
(QMC),— the T- matrix approximation ) -'^^i -'^^ a general- 
ized RPA)2^ the parquet approachj^i the temperature- 
cutoff renormalization-group (TCRG)r^ and the two- 
particle self-consistent (TPSC) approximation.^ Most of 
these approaches are limited to the paramagnetic state 
and (or) at the van-Hove filling. 

It is therefore of interest to theoretically investigate 
magnetic fluctuations in a low-dimensional band ferro- 
magnet by incorporating correlation effects within an 
approach in which spin-rotational symmetry and the 
Goldstone mode are explicitly preserved. The inverse- 
degeneracy expansion scheme^'* extended recently to 
study band ferromagnetism in the Hubbard mode l^^i^^ 
is particularly suitable. In this approach correlation ef- 
fects are systematically incorporated by including self- 
energy and vertex corrections so that spin-rotation sym- 
metry and Goldstone mode are preserved order by or- 
der. It therefore provides a quantitative description of 
the long- wavelength, low-energy magnetic fluctuations 
consistent with the continuous spin-rotational symmetry, 
which play a dominant role in determining the magnetic 
behaviour of low-dimensional systems such as ultrathin 
films and magnetic nanostructures. 

Correlation effects incorporated beyond RPA were 
shown to strongly suppress ferromagnetism, as quan- 
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titatively demonstrated for several three-dimensional 
latticesj2^ Furthermore, the correlation-induced spin- 
charge coupling was shown to provide a strong magnon- 
damping mechanism for modes lying within the Stoner 
gapf^^ whereas magnon damping within RPA is only 
due to decay of collective excitations into single parti- 
cle Stoner excitations. 

In this paper we will quantitatively investigate band 
ferromagnetism within the t-t' Hubbard model on a 
square lattice. We will evaluate quantum corrections 
to the magnon propagator by incorporating correla- 
tion effects in the form of self-energy and vertex cor- 
rections within the Goldstone-mode preserving inverse- 
degeneracy expansion scheme. This will allow us to ex- 
amine the stability of the ferromagnetic state with re- 
spect to both long- and short-wavelength modes. We 
will show that the enhanced ferromagnetic stability near 
the van Hove filling for large t' (~ 0.5) can be read- 
ily understood in terms of strongly suppressed saddle- 
point contribution to the exchange energy term in the 
spin stiffness. Finally, we will investigate the finite tem- 
perature electron spin dynamics directly in terms of the 
electronic spectral-weight transfer across the Fermi en- 
ergy, and thereby determine the thermal magnetization 
decay. In order to highlight the Goldstone-mode preserv- 
ing character of our approach which is in accord with 
the Mermin- Wagner theorem, we will focus on the con- 
tribution of long-wavelength modes and determine the 
thermal magnetization decay and Tc in the presence of a 
small anisotropy gap. 

The organization of this paper is as follows. In Sec. 
II we review the Goldstone-mode preserving diagram- 
matic approach for quantum corrections to the trans- 
verse spin fluctuation propagator. Sections III and IV 
discuss the stability of ferromagnetic state with respect 
to long- and short-wavelength fluctuations in terms of 
the spin stiffness and a characteristic magnon energy for 
zone-boundary modes, respectively. In Section V we dis- 
cuss finite-temperature effects in terms of spectral-weight 
transfer across the Fermi energy due to electron-magnon 
coupling, and the resulting magnetization decay and Tc 
in the presence of a small anisotropy gap. Finally, con- 
clusions are presented in Section VI. 



II. TRANSVERSE SPIN FLUCTUATIONS 

Transverse spin fluctuations are gapless, low-energy 
excitations in the spontaneously-symmetry-broken state 
of magnetic systems possessing continuous spin-rotation 
symmetry, and play an important role in diverse macro- 
scopic properties such as existence of long-range magnetic 
order, temperature dependence of magnetization, tran- 
sition temperature, spin correlations etc. We therefore 
consider the time-ordered, transverse spin-fluctuation 
propagator 
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FIG. 1: First-order quantum corrections to the irreducible 
particle-hole propagator. 



in terms of the electron spin-lowering and -raising opera- 
tors — 5'|(cr^/2)'I'i, evaluated in the spontaneously- 
symmetry-broken ground state |5'g)- For simplicity, we 
consider the saturated ferromagnetic state, in which the 
minority-spin states lie above the Fermi energy; the mag- 
netization m is hence equal to the particle density n. 

The spin-fluctuation propagator in q, a; space can be 
expressed generally as 



l-C/0(q,c^) 



(2) 



in terms of the exact irreducible particle-hole propagator 
(/)(q, uj), which incorporates all self-energy and vertex cor- 
rections. The inverse-degeneracy (l/Af) expansion2ii2^ 

^ = (3) 

systematizes the diagrams for (/) in powers of the expan- 
sion parameter 1 /Af which, in analogy with 1/ S for quan- 
tum spin systems, plays the role of h. Required from 
spin-rotation symmetry, the exact cancellation (p^P^^ = 
of all quantum corrections {p > 1) for q,uj — en- 
sures that the Goldstone mode is preserved order by or- 
der. From (2) it follows that U^(j)''P\q, Lu) represents the 
magnon self energy at order p. 

The bare particle-hole propagator is given by 



X^{t~n^^{^o\T[Sr{t)S+{t') 



(1) 



</>(o)(q,c.)^x''(q,c^) = E' 



k-q 



11] 



(4) 
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in terms of the HF-level ferromagnetic band energies 
= ek — crA and the exchange band spHtting 2 A = mil; 
the superscript +(— ) refers to particle (hole) states above 
(below) the Fermi energy ep. As only the "classical" term 
0(0) survives in the Af ^ oo limit, the RPA ladder sum 
X°(q, cj)/l — C/x°(q, cj) amounts to a classical (unrenor- 
malized) description of non-interacting spin-fluctuation 
modes. 



The first-order quantum corrections 4>^^\ involving 
self-energy and vertex corrections of order 1/Af, have 
been obtained recently for a saturated ferromagnet,— 
and physically incorporate such effects as minority-spin 
spectral-weight transfer, correlation-induced exchange 
correction, and coupling of spin and charge fluctuations. 
Shown diagrammatically in Fig. 1, the quantum correc- 
tions were obtained as: 
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The exact cancellation of for q,^ — 0, required 
from spin-rotation symmetry in order to preserve the 
Goldstone mode has been discussed earlier. Indeed, the 
cancellation holds for all ui, indicating no spin- wave am- 
plitude renormalization, as expected for the saturated 
ferromagnet in which there are no quantum corrections 
to magnetization. 

For flnite q, the magnon energy Wq is obtained from the 
pole condition 1 — Udi<j){<l, — Wq) = in Eq. (2). Evalua- 
tion of the f2-integral in (5) was carried out numerically 
by including the contribution of both the low-energy col- 
lective spin-wave excitations as well as the high-energy 
Stoner excitations, as discussed earlier in detail. 

We consider the Hubbard model on a square lattice 
with band dispersion 

eic = — 2i(cos fcj; + cos ky) + At' cos cos ky (6) 

corresponding to NN and NNN hoppings t and t' . It is 
instructive to expand the dispersion (6) around special 
points in k space. Near (±7r,0), the dispersion exhibits 
the saddle-point behaviour 

e],^{t + 2t')nl-{t-2t')nl-At' (7) 



in terms of the shifted momentum k, with a similar be- 
haviour near (0, ±7r) with and Ky interchanged. This 
saddle-point behaviour results in the characteristic loga- 
rithmic density-of-states (DOS) singularity 

7V(e)^ln^— (8) 

at energy —At' , as seen in Fig. 2. In the limit t' 
t/2, the effectively one-dimensional dispersion results in 
a stronger (~ + 2t) DOS singularity at the band 

bottom. Near the origin (0, 0) we have 

ek = (i - 2t')\^ + t'klkl - A{t - t') , (9) 

and the quartic dispersion surviving in the limit 
t' t/2 yields a marginally stronger DOS singularity 
(l/Ve + 2t) In at the band bottom. 

We choose an energy scale such that t = 1. The effects 
of the distribution of spectral weight in DOS have been 
investigated for < <' < 0.5 in the intermediate range 
of interaction strengths {U/W = 0.5, 1.0 and 1.5), where 
the electron bandwidth W = 8t is independent of t'. 
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FIG. 2; DOS for band dispersion (6), showing the enhanced 
saddle-point contribution and van Hove singularity with in- 
creasing next-nearest-neighbor hopping t' . 



III. SPIN STIFFNESS 

The spin stiffness D = ojq/q^ in the ferromagnetic 
state, defined in terms of the magnon energy LUq for small 
q, provides a quantitative measure of the stability of 
the ferromagnetic state against long-wavelength fluctua- 
tions, with negative D signalling loss of long-range mag- 
netic order. We first review the different contributions to 
spin stiffness at the RPA (classical) level as their behav- 
ior and interplay provide insight into the magnitude of 
the first-order quantum corrections and of ferromagnetic- 
state stability, as discussed in the next subsection. 



A. RPA level: hopping and exchange contributions 

Expanding x''(q, lj) for small q, w the RPA (classical) 
spin stiffness can be expressed as^ 



d 



2A 



(10) 



in d dimensions. Here the angular bracket ( ) represents 
momentum summation normalized over the number of 
occupied states (;^X]k)- '^^^ terms in (10) of or- 
der t and t'^/U represent hopping and exchange contri- 
butions to the spin stiffness, respectively, corresponding 
to delocalization-energy loss and exchange-energy gain 
upon spin twisting. 

The stability of the ferromagnetic state therefore in- 
volves a subtle competition between the hopping con- 
tribution which favors the FM state, and the exchange 
contribution which tends to destabilize it. Both the hop- 
ping and exchange contributions are sensitive to band 
dispersion and filling, as discussed below. 
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FIG. 3: Hopping (upper left) and exchange (upper right) con- 
tributions to the RPA spin stiffness (lower) at U/W — 1.5, 
showing enhancement in stiffness and reduction in optimal 
density with increasing t' . 



Figure 3 shows the hopping and exchange contribu- 
tions to the spin stiffness. The limit n ^ 1 of filled 
majority-spin band is physically simple and essentially 
independent of band-dispersion. While the hopping con- 
tribution to spin stiffness vanishes, the exchange contri- 
bution approaches a finite value, thereby destabilizing 
the ferromagnetic state. 

The low-density (n <C 1) behavior is, however, highly 
sensitive to the band dispersion and the distribution of 
spectral weight near the band bottom. With increas- 
ing t' , although the hopping contribution decreases and 
nearly saturates, the exchange contribution decreases 
rapidly and even vanishes as t' approaches 0.5, result- 
ing in substantial enhancement in the stability of the 
ferromagnetic state. This behaviour of the hopping and 
exchange contributions can be readily understood from 
key features of the band dispersion, as explained below. 

For low filling, the occupied states he near k — (0, 0), 
so that from Eq. (6) 



V2 



cos ky) — 8t' cos kx cos ky 



2i(cos kx 
~ (4i - 8t') ~{t~ W)k^ , 

leading to a delocalization contribution 
1 



(11) 
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[V^e^)^{t-2t')-{t-U')2iTn/A (12) 



to spin stiffness which decreases linearly with filling n. 
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FIG. 4: The corrected spin stiffness shows a strong enhance- 
ment in the iow density regime as t' approaches 0.5, with 
optimai density approaching the van Hove fiiiing. 



with the slope changing sign at t' = t/4, as seen in Fig. 
3. Here we have used n = kp/Air corresponding to the 
circular filling pocket. 

Similarly, for the exchange term in the low-density 
limit, we obtain 

(Vek)^ = (2t - At' cos kyf sin^ + x ^ y (13) 
{2t - 4t')^[k2 - {kt + fc^)/3] + 8t'{2t - At')klkl , 

which yields an exchange contribution 

27rri / 2t' 



(0) _ 47rft - 2t' 

cxch jj 



1 - 



1 - 



t-2t' 



(14) 



decreasing quadratically with [t — 2t') and linearly with 
n, with the slope changing sign at t' = t/4. 



Close to van Hove filling, what is the contribution of 
states near the DOS singularity to the delocalization and 
exchange contributions to spin stiffness? As the DOS sin- 
gularity is associated with saddle points where Vek = 0, 
states near the saddle points necessarily yield a small 
exchange contribution and a finite delocalization contri- 
bution. With increasing t', the fraction of states near 
the saddle points increases, and the exchange contribu- 
tion correspondingly decreases rapidly, as seen in Fig. 3. 
Indeed, for t' = t/2, the gradient Ve^ vanishes for all 
states on the two axes kr^ — and fcj, = 0, as seen from 
(13), and as all these states lie at the band bottom {~2t), 
the exchange contribution vanishes as n — > 0. 

With increasing U, both the density range and stiff- 
ness of FM state increase due to the smaller exchange 
contribution, while the hopping contribution and hence 
the optimal density remains unaffected. Similar effects 
of U and t' on ferromagnetic tendency have been found 
by considering the instability of the paramagnetic state 
within the generalized RPAi^ 

As explicitly shown earlier j^^i^^ the first-order quan- 
tum corrections to spin stiffness essentially involve 
correlation-induced exchange processes, and therefore 
the behavior of classical exchange contribution actually 
provides a qualitative idea of the first-order quantum cor- 
rections as well. The above analysis therefore provides 
insight into the origin of the enhancement of ferromag- 
netism near the van Hove filling at sufhciently large t' . 



B. Quantum correction 

The expansion of the first-order quantum correction 
(f)^^") for small q has been carried out earlier j^^i^^ and 
yields the corresponding quantum correction to spin stiff- 
ness 
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All four terms in (15) represent correlation-induced ex- spectral-weight transfer is a correlation effect corre- 
change processes involving minority-spin intermediate sponding to the possibility of a site being unoccupied 
states which are transferred down in energy. The by a majority-spin electron. Quantum correction to 
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FIG. 5: The stiffness as well as the stable range of densities 
are enhanced with increasing U while optimal density remains 
unaffected. The ferromagnetic state is unstable for t' < 0.35 
at all considered values of U, for t' = 0.4 and U/W < 1.0, 
and also for t' = 0.45 at U/W = 0.5. 



FIG. 6: Characterizing large-momentum, zone-boundary 
modes, the renormalized (ren) dominant-mode energy shows 
a strong suppression at low densities due to quantum correc- 
tions, restricting stability to an intermediate range of densi- 
ties, shown for U/W = 1.0. 



the delocalization contribution was shown to vanish 
identically j^^'^^ Together with the classical exchange con- 
tribution (10), the above correlation-induced exchange 
contribution yields a further reduction in the renormal- 
ized spin stiffness D = - D^^). 

Overall, we find that quantum corrections are sup- 
pressed by t\ as expected from the decreasing classi- 
cal exchange contribution (Fig. 3). Nonetheless, quan- 
tum corrections are strong enough to destabilize the FM 
state for t' < 0.35 for all considered interaction strengths 
and in the whole range of electronic density. This is in 
excellent quantitative agreement with results of the T- 
matrixiS and generalized-RPA-" approaches. 

Figure 4 shows the D vs. n behaviour. There is a 
significant enhancement of ferromagnetism in the low- 
density regime as t' approaches 0.5. While optimal den- 
sity decreases and approaches van Hove filling, the sta- 
ble density range increases due to significant enhance- 
ment of spin stiffness in the low density regime. This en- 
hancement of ferromagnetism with increasing t' is in gen- 
eral agreement with results obtained from various other 
approachesM^^^-^^ 



Figure 5 highlights the interplay of t' and U. With 
increasing t' , the ferromagnetic state becomes stable even 
at relatively small interaction strength. The FM state 
is unstable in the whole density range for t' = 0.4 at 
U/W < 1.0 and similarly for t' = 0.45 at U/W = 0.5, but 
becomes stable in a limited density range at intermediate 
coupling [U/W— 1.0 and 1.5), with an optimization near 
n — 0.3. For t' — 0.5, quantum correction becomes less 
significant, resulting in stabilization of the ferromagnetic 
state even at U/W = 0.5, and significant enhancement 
in spin stiffness with decreasing n. 

For both t' = 0.45 and 0.5, increasing U results in en- 
hancement of spin stiffness as well as of the stable range of 
densities. This is because of both enhanced RPA stiffness 
due to reduced exchange contribution and also reduced 
quantum corrections. 

Our results for the stability of the ferromagnetic 
state for t' — 0.4, 0.45, and 0.5 are in good quanti- 
tative agreement with the prediction of the T-matrix 
approximatior>i^ in the low-density regime which, how- 
ever, predicts a relatively large upper density limit. This 
quantitative difference with respect to the upper density 
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limit may be attributed to the over- estimation by the T- 
matrix approach which is more reHable in the low density 
limit.— Furthermore, our results for t' — 0.5 are in 
excellent quantitative agreement with the predictions of 
the generalized RPA.^° Our overall finding of decreasing 
critical value of t' with increasing \J beyond which the 
ferromagnetic state is stable is in good qualitative agree- 
ment with the approaches mentioned in the Introduction. 



IV. MAGNON-ENERGY RENORMALIZATION 
FOR LARGE-MOMENTUM MODES 
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FIG. 7: The first-order correction to spin-^ electron Green's 
function due to electron-magnon interaction, which results in 
spectral-weight transfer across the Fermi energy and thermal 
magnetization decay. 



We now consider quantum corrections for zone- 
boundary modes with energy Wq near the magnon DOS 
peak energy (w*). The behavior of this "dominant-mode" 
energy with electron filling n is of interest as it deter- 
mines the stability of the FM state with respect to typical 
short- wavelength fluctuations. We have checked for sev- 
eral such q and find the behavior to be quantitatively 
very similar. 

Figure 6 shows the behavior of the renormalized 
dominant-mode energy with electron filling for t' = 0.45 
and 0.5, along with the RPA results. In contrast to Fig. 
5, the nearly similar behavior shows the relatively much 
weaker sensitivity of large-q modes on details of band 
dispersion. The strong suppression at low densities, a 
trend already present at the RPA level, results in an 
intermediate-density range of stability with respect to 
short-wavelength fluctuations. 

Comparison of Figs. 5 and 6 shows a marked contrast 
in the behavior of magnon energy for long- and short- 
wavelength modes, especially for t' — 0.5. While positive 
spin stiffness in the low-density range indicates stability 
with respect to long- wavelength fluctuations, the van- 
ishing of dominant-mode energy w* for relatively short- 
wavelength modes signals the spontaneous onset of strong 
local fluctuations. On the other hand, while to* shows a 
peak near n « 0.6, indicating stability with respect to lo- 
cal fluctuations, the vanishing spin stiffness near n — 0.6 
indicates loss of long-range magnetic order. That the sta- 
bility of the ferromagnetic state with respect to long- and 
short-wavelength fluctuations can differ so substantially 
is a significant feature of our approach. 

The substantial decrease in magnon energy at large 
momenta in the low-density regime (Fig. 6) is of 
considerable relevance in the context of experimen- 
tally observed "anomalous softening" in ferromagnetic 
manganites.^- 

In the intermediate-density range where FM state 
is locally stable, but unstable with respect to long- 
wavelength fluctuations, some commensurate spin den- 
sity wave states other than ferromagnetic and antiferro- 
magnetic were found in the generalized RPA2^ investiga- 
tion. 



V. SPECTRAL- WEIGHT TRANSFER AND 
THERMAL DECAY OF MAGNETIZATION 

As spin-rotation symmetry and the Goldstone mode 
are explicitly preserved in our approach, it is therefore in 
accord with the Mermin- Wagner theorem. Long wave- 
length modes in one and two dimensions yield divergent 
thermal fluctuations, indicating absence of long-range 
magnetic ordering at any finite temperature. With quan- 
tum corrections included to yield the correct magnon 
energy scale for long wavelength modes, our approach 
therefore provides a quantitative tool for investigating 
fluctuation effects in low-dimensional band ferromagnetic 
structures. As a simple illustration, we consider in this 
section the reduction in magnetization with temperature 
due to thermal excitation of magnons, suppression of 
the infrared divergence by anisotropy induced gap in the 
magnon spectrum, and the consequent stabilization of 
the ferromagnetic state in two dimensions. 

We first consider a simple description of the magnon 
propagator at finite temperature. Including both quan- 
tum and thermal corrections, the magnon propagator is 
approximately given by 



x-+(Q,f^) 



(16) 



where 



?TIQ = mc^.{Sz)T/{Sz)o 
= nci.{Sz)T/{Sz)o 



(17) 



are the renormalized magnon amplitude and energy, with 
()t denoting the thermal average at temperature T. 
These thermal renormalizations arise from the 0{lu) term 
in the quantum corrections, which is proportional to the 
magnetization {n^)T ^ ("i)T = '2{Sz)t- At zero tem- 
perature, quantum corrections to the 0{lo) term cancel 
identically, as expected from the absence of any quan- 
tum correction to magnetization. The above form of 
the magnon propagator, similar to that for a quantum 
Heisenberg ferromagnet, ensures that the sum-rule 



{[S+,S-])t = 2(Sz)i 



(18) 
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following from the spin commutation property is identi- 
cally satisfied at all temperatures. The gapped Stoner 
excitations, which distinguish a band ferromagnet, have 
been neglected in the above consideration as our focus is 
on long- wavelength, low-energy fluctuations. 

The thermal reduction of magnetization in a band fer- 
romagnet arises from the electron spectral-weight trans- 
fer across the Fermi energy. Due to the spin-flip scatter- 
ing of electrons accompanying the thermal magnon exci- 
tations, a portion of the spin-l spectral weight is trans- 
ferred to the spin-| band above the Fermi energy, while 
an equal amount of spin-| spectral weight is transferred 



to the spin-l band below the Fermi energy. The corre- 
sponding changes in the electron densities due to the first 
order fluctuation process (Fig. 7) are obtained as: 



5n^ = -5n^ = / — Iin[G;(k, c^)I]T(k, c^)Gj(k, c^)] , 

(19) 

where GJ is the HF-level (advanced) propagator and the 
electron self-energy correction due to electron-magnon 
coupling is given by 



sT(k,a;) = U^Yl I |^X-+(Q,f^)Gi(k~Q,u;-17)[- 



-1 



Q 



(20) 



Evaluating the spectral weight transfer from Eq. (19), 
we obtain the thermal reduction in magnetization 

{Sz)o - {Sz)t = (Sni - Sn^)/2 = Sn^ 

^ u^y-^^y-f— — i — _V (21) 

Coupled with Eq. (17) for magnon amplitude and energy 
renormalization, the above equation self-consistently de- 
termines the magnetization {Sz)t at temperature T. The 
dependence of the magnon energy naturally yields 
divergent reduction in one and two dimensions, and the 
low-temperature Bloch-law (T^^^) fall off in three dimen- 
sions. 

As the temperature approaches Tc from below, both 
the magnetization {Sz)t and the renormalized magnon 
energy CIq become vanishingly small, and from Eq. (21) 
we therefore obtain 




where we have taken toq w 2{Sz)o and ei(._Q — e]^ ~ 
2A = nil for the dominant small-Q modes. With the 
formal identification n = 2S, this expression is similar 
to that obtained by mapping the itinerant electron sys- 
tem to an effective Heisenberg model.— However, to our 
knowledge, this is the first such result for an itinerant fer- 
romagnet which has been obtained directly in terms of 
the spectral weight transfer. 



For a two-dimensional isotropic system with Qq ^ 
Q^, the transition temperature vanishes due to the log- 
arithmically divergent contribution of long-wavelength 
magnon modes. Presence of an anisotropy gap in the 
long-wavelength magnon spectrum flq = DQ^ + Aq due 
to intrinsic spin anisotropy stabilizes the ferromagnetic 
ordering at finite temperature by suppressing the di- 
vergent contribution of long-wavelength modes. Recent 
magnetization measurements^ on ultrathin Fe films have 
indeed found substantial enhancement in thermal exci- 
tation of magnons with decreasing anisotropy gap which 
is sensitive to the interface between film and substrate. 
With an anisotropy gap, the expression for transition 
temperature reduces to 

T,^nD/\ni^^j , (23) 

where A ^ 1 is an upper-momentum cutoff. This re- 
sult for anisotropy-driven stabilization of ferromagnetic 
ordering at finite temperature is in close analogy with 
similar investigations for localized spin systems.— With 
quantum corrections included in the spin stiffness D, the 
above result highlights the role of correlation effects on 
long-wavelength fluctuations in essentially determining 
the transition temperature. 

VI. CONCLUSIONS 

In conclusion, we have investigated quantum cor- 
rections to collective magnon excitations in a two- 
dimensional correlated band ferromagnet. Self-energy 
and vertex corrections were incorporated within a spin- 
rotationally-symmetric scheme in which the Goldstone 
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mode is explicitly preserved. By studying the ferromag- 
netic state spin stiffness, we showed that ferromagnetism 
is strongly enhanced near the van Hove filling for t' ~ 0.5, 
in agreement with earlier studies from the paramagnetic 
side. This enhancement in ferromagnetism was shown 
to be a consequence of strongly suppressed saddle-point 
contribution near the van Hove singularity to the com- 
peting exchange part of the spin stiffness. 

Due to the strong spin-charge coupling, resulting in 
an anomalous momentum dependence of the magnon 
self energy, magnon energies for zone-boundary modes 
were found to be strongly suppressed in the low-density 
regime. Also obtained earlier in our three-dimensional 
calculations^^ this magnon-energy suppression for zone- 
boundary modes is relevant in the context of experimen- 
tally observed "anomalous softening" in ferromagnetic 
manganites.'^ 

We also studied the finite-temperature electron spin 
dynamics due to coupling with the fluctuating trans- 
verse field of collective excitations in the presence of a 
small anisotropy gap. The magnetization decay and Tc 
were obtained directly in terms of the electronic spec- 
tral weight transfer, rather than first mapping the itiner- 



ant ferromagnet to an effective Heisenberg spin model. 
The renormalized spin stiffness was shown to essen- 
tially determine the energy scale for magnetization de- 
cay and Tc- The vanishing of Tc in the isotropic limit 
due to divergent long wavelength fluctuations is in ac- 
cord with the Mermin- Wagner theorem. The Goldstone- 
mode preserving approach thus naturally allows for long- 
wavelength, low-temperature, and low-dimensional stud- 
ies in the broken-symmetry state. 

These results are relevant in the context of recent mag- 
netization measurements on ultrathin transition-metal 
films, although quantitative comparison with experi- 
ments would necessitate incorporation of several realistic 
features such as multiple 3d orbitals as well as presence 
of multilayers and non-magnetic substrates. 
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